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Abstract 

We present optical linear spectropolarimetry of the microquasar LS I +61° 303. 



The continuum emission is mildly polarized (up to 1.3 %) and shows almost no tem- 
poral change. We find a distinct change of polarization across the Ha emission line, 
indicating the existence of polarization component intrinsic to the microquasar. We 
estimate the interstellar polarization (ISP) component from polarization of the Ha 
line and derive the intrinsic polarization component. The wavelength dependence of 
the intrinsic component is well explained by Thomson scattering in equatorial disk of 
the Be-type mass donor. The position angle (PA) of the intrinsic polarization ~ 25° 
represents the rotational axis of the Be disk. This PA is nearly perpendicular to the 
PA of the radio jet found during quiescent phases. Assuming an orthogonal disk-jet 
geometry around the compact star, the rotational axis of the accretion disk is almost 
perpendicular to that of the Be disk. Moreover, according to the orbital parameters of 
the microquasar, the compact star is likely to get across the Be disk around their pe- 
riastron passage. We discuss the peculiar circumstellar structure of this microquasar 
inferred from our observation and possible connection with its high-energy activities. 

1. Introduction 

A microquasar is an X-ray binary system displaying relativistic radio jets (see review 
of Ribo 2005 and references therein). It consists of a compact star (neutron star or stellar- 
mass blackhole) and a mass donor star. About 20 microquasars have been identified in our 
galaxy to date (Paredes 2005). They are likely to be scaled-down version of radio- loud active 
galactic nuclei because of their morphological and physical similarities (relativistic jet ejection 
and blackhole engine). A microquasar provides an excellent laboratory for study of mass 
accretion and ejection phenomena in a strong gravitational field on human timescales. Their 
spectra extend from radio to 7-rays (e.g. Mirabel & Rodriguez 1999). Optical emission from 
microquasars are usually considered to be either thermal emission emerged from outer part of 
the accretion disk around the compact star (low-mass donor case) or stellar light of the mass 
donor (high- mass case). However, recent observations suggest that some microquasars show 
synchrotron emission even at optical wavelengths in outburst phases (Kanbach et al. 2001; 
Uemura et al. 2004). 

Optical polarimetric observation can be unique to probe microquasar systems, since 
we can obtain information on jets (e.g., geometry, magnetic field, electron energy) from syn- 
chrotron emission or on geometrical properties of the binary system from the light scattered 
by circumstellar matter (accretion disk and/or stellar wind) in a scale not being accessible by 
contemporary imaging techniques. Several optical polarimetric studies of microquasars have 
been reported for, e.g., LS 5039 (Combi et al. 2004), GRO J1655-40 (Scahriti et al. 1997; 
Gliozzi et al. 1998), and SS 433 (Efimov et al. 1984) so far. They revealed the presence 
of intrinsic polarization. However, all of these studies employed broad-band filters which can 
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not reveal how the polarization correlates with the spectral features. Spectropolarimetry per- 
mits separation of the line and continuum components, providing a more accurate way to find 
intrinsic and interstellar components of the observed polarization. 

We made optical spectropolarimetric observations of the microquasar LS I +61° 303. 
This is one of Be/X-ray binary systems, having a large eccentricity of e = 0.72 ±0.15 (Casares 
et al. 2005). This object is an optical counterpart of the variable radio source GT 0236+610, 
which was discovered within la error circle of the 7-ray source 2CG 135+01 = 3EG J0241+6103 
(Gregory & Taylor 1978). This source exhibits radio outbursts every 26.5 days (Taylor & 
Gregory 1982). Nearly identical periodicity has been found in X-rays (26.42 ± 0.05 days; 
Paredes et al. 1997). They are consistent with the orbital period determined by optical and 
infrared observations (Hutchings & Crampton 1981; Paredes et al. 1994; Gregory 2002). 
Similar periodicity has also been suggested in 7-rays (27.4 ± 7.2 days; Massi 2004b), which 
supports the physical connection of LS I +61° 303 with 3EG J0241+6103. Recently, Albert 
et al. (2006) found very high energy 7-ray emission (above 100 GeV) showing a possible 
periodicity. Based on the lUE ultraviolet and ground-based spectroscopy, the mass donor of 
LS I +61° 303 is classified as a BOVe star (Hutchings & Crampton 1981). Long monitoring 
of radio flux (Gregory et al. 1999) and Ha line emissions (Zamanov et al. 1999) indicated 
that the mass-loss rate of the Be star has a period of about 4.6 yr. Radio interferometric 
observations found that the axis of radio jet rapidly rotated in the projected sky plane (Massi 
et al. 2004a), giving an S-shaped morphology like SS 433 (Hjellming & Johnston 1988) and a 
'micro-blazer' model (Kaufman Bernado, Romero, Mirabel 2002; Romero, Kaufman Bernado, 
& Mirabel 2002) is proposed for this binary. 

In this paper, we examine a circumstellar structure of LS I +61° 303 with our new 
spectropolarimetric data. We describe the observations and data reduction in §2, and show 
the results in §3. We discuss the circumstellar structure of LS I +61° 303 in §4. Then, we 
summarize our study in §5. 

2. Observations and Data Reduction 

The observations were performed on two nights in 2005 January (epoch 1) and on six 
nights in 2005 October and November (epoch 2) with a low-resolution spectropolarimeter, 
HBS (Kawabata et al. 1999), attached to the 1.88m telescope at Okayama Astrophysical 
Observatory, National Astronomical Observatory of Japan (NAOJ). The log of the observations 
is shown in Table 1. HBS has a superachromatic half- wave plate and a quartz WoUaston prism, 
giving a capability of measuring linear polarization in a wide range of optical region (4000-9000 
A). We employed a SITe 512 x 512 x 24 CCD camera (MicroLuminetics, Inc.) with in 
epoch 1 and a Marconi 2048 x 512 x 13.5 fim CCD camera (Andor Tech, Inc.) in epoch 2. We 
used a slit of 0.2mm width, yielding a wavelength resolution of ~ 60 A. A unit of the observing 
sequence consists of successive integrations at 0°, 22.5°, 45° and 67.5° position angles of the 
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half- wave plate. 

Data reduction was performed using a standard software developed for HBS observations, 
HBSRED (Kawabata et al. 1999). For polarimetric calibration, we used data of polarized and 
unpolarized standard stars (Wolff, Nordsieck, & Nook 1996), obtained with or without a fully- 
polarizing, Glan- Taylor prism. Polarimetric accuracy estimated from peak-to-peak variation of 
unpolarized standard star data was ~ 0.05% over the observed wavelength range. The flux was 
calibrated using data of spectrophotometric standard stars (Taylor 1984). 

3. Results 

3.1. Overall polarization properties 

Figure 1 shows observed polarization spectra (degree of polarization Pobs and its position 
angle 6'obs) on 2005 Oct 27 and Nov 8. The continuum emission exhibits rather smooth wave- 
length dependence over the optical region. The nightly-averaged V-band polarization data are 
shown in Tablel, which indicates that there was no significant time variation in the continuum 
polarization over the period of our observations. Figure 2 shows an enlarged plot around the 
Ha emission line. It indicates an apparent change of polarization across the Ha line, suggesting 
the existence of a polarization component intrinsic to LS I -|-61° 303. We discuss it in §3.3. 

3.2. Equivalent width of Ha emission line 

Ha emission of a Be star is thought to originate mostly from equatorial disk around the 
star (e.g., Quirrenbach et al. 1997). Paredes et al. (1994) and Zamanov et al. (1999) revealed 
that equivalent width (EW) of Ha line in LS I -|-61° 303 changes with a possible period of ~4.6 
yr. This is considered to reflect the variation of mass-loss rate of the Be star. Figure 3 displays 
the mass-loss phase dependency of the EWs derived from our data together with those from 
literature, which shows that EWs of our data were commonly at low level compared with the 
past observations. It suggests that the disk of the Be star has not developed very much and 
the ~4.6 yr periodicity is unclear during our observation. According to AAVSO^ database any 
significant brightness change was not reported during the period of our observation. These 
facts imply that the binary system was in a quiescent phase in the optical light. 

3.3. Interstellar polarization 

Observed polarization is generally expressed as a vectorial sum of intrinsic polarization 
and interstellar polarization (ISP), i.e., 

gobs (A) = (?int(A) + Pisp(A) cos26'isp, Uobs(A) = Mint (A) + Pisp(A) sin26'igp, (1) 

where (gobs, ^obs) and (gint, ^^int) are Stokes parameters of the observed polarization and the 
intrinsic polarization, respectively. Pjsp and 9isp are degree of the ISP and its position angle, 

^ http://www.aavso.org/ 
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respectively. ISP is produced by dichroic absorption by magnetically-aligned aspherical dust 
grains existing between the object and the earth. To derive the intrinsic polarization, it is 
essentially important to know the ISP component accurately. The degree of ISP is well expressed 
by a smooth function of wavelength in optical region as, 

Pisp(A) = P^axexp(-i^ In^^), (2) 

where Pmax is the peak polarization degree of the ISP at wavelength Amax (Serkowski et al. 
1975). Upper-limit of Pmax is practically fixed by color excess, E^.y, as 

Pmax < 9.0 X Eb^v (3) 

(Serkowski et al. 1975). Paredes & Figueras (1986) derived Eb-v = 1-13 for LS I +61° 303, 
which leads to Pmax < 10.2% and a median value of Pmax ~ 5.1% (Serkowski et al. 1975). The 
observed polarization degree Pobs.y ~ 1-3% (See Table 1) is much smaller than the expected 
degree of polarization. It suggests that the ISP is partly canceled out by non-zero intrinsic 
polarization of which the PA is nearly orthogonal to that of the ISP. To confirm it, we fit the 
Serkowski function (Eqn. 2) to the observed polarization data and derived the ISP parameters. 
Figure 4 shows the obtained K and Amax values, together with the i^-Amax relation derived by 
Whittet et al. (1992), 

K = (1.66 ± 0.09)Amax(/im) + 0.01 ± 0.05 (4) 

This figure shows that most of our data points detach from the linear /C-Amax relation, which 
is also consistent with the idea that the observed polarization has a substantial intrinsic com- 
ponent. 

To evaluate the ISP component, we used the change of polarization across the Ha emis- 
sion line. Line emission generally emerges from optically-thin regions and is subject to less 
scattering than continuum emission. This leads to a lower polarization level in the line emis- 
sion than in the continuum emission. In the case of LS I +61° 303, the line emission is thought 
to arise from ionized gas of the Be disk. Without a few exceptions, we can consider that the 
line emission in Be stars should show little or no polarization (e.g. Poeckert 1975; Poeckert 
& Marlborough 1977; McLean & Clarke 1979; Quirrenbach et al. 1997). Therefore, we 
assume that the polarization of the line component practically represents the ISP. To subtract 
the continuum component from the polarization at Ha emission line (denoted by suffix 'em'), 
we used its neighboring continuum polarization on both sides of the line ('conl' and 'con2'). 
We derive Stokes parameters I,Q,U of the line emission as 

/ AA ( ^'^^ -^conl ~l~ -^con2 \ 

vAAem AAconl + AAcon2/ 

Q / Qem ^ -^em 'Zconl ^ -^conl + 900112 ^ -^con2 \ 

V AAcm AAconl + AAcon2 / 

jY A \ f ^cm ^ -^cm ^conl ^ -^conl ~l~ ^con2 ^ -^con2 \ /iy\ 

Uua = AAem I — TT ^ — TT I , (7j 

\ AiAcm ^^conl + ^^con2 / 



where AA is the band width of each component (em, conl and con2). We set all AA to be 150 



component in those bands except for the Ha line. We calculate the polarization degree and PA 
as 



The derived Ha polarization is shown in Tablel. Weighted-mean and standard deviation of 
all the eight nights data are PHa=2.16 ± 0.20% and ^Ha=126.4° ± 4.5°, respectively. These 
values are comparable to the polarization of the field stars within 3 degree (4.3 ± 1.1 %, 112 
± 7°; Heiles 2000). Since PA of the ISP is almost constant with wavelength, we consider 
that represents ^isp. To derive the spectrum of the intrinsically polarized component, we 
should know the wavelength dependence of degree of the ISP, further. Once Amax is fixed, the 
remaining ISP parameters, K and Pmax? are calculated from Equations (2), (4) and the derived 
Pho- Serkowski et al. (1975) indicated that ISP has a peak polarization at ~ 5400 ± 600 
Ain the case of Pmax/ Eb-v < 3. We thus assume that Amax is 5000 A, 5500 A, or 6000 A, and 
calculate the intrinsic polarization. 

3.4- Origin of intrinsic polarization 

Figure 5 shows the spectra of the intrinsic polarization for different Amax- The wavelength 
dependence of the intrinsic polarization is closely similar to those of typical Be stars, e.g. ( 
Tau (Wood et al. 1997). The polarization of Be stars is explained by the result of Thomson 
scattering in their equatorial disks. If there is only one polarizing component (e.g. a single disk), 
the PA of the intrinsic polarization 6[j^t should be constant against wavelength. This assumption 
is reasonable for most Be stars. The derived 6'int becomes rather flat for Amax, suggesting that 
the assumptions of the range of Amax and the Be star-origin polarization are self-consistent. 
Therefore, we conclude that the origin of the intrinsic polarization is predominantly Thomson 
scattering in the Be star. 

4. Discussion 

4-1- Topology of the Be equatorial disk and orbital plane 

Neither 6'int nor the EW of the Ha line do show temporal variation, which suggests that 
the Be equatorial disk is nearly stable during the period of our observation, di^t ~ 25° implies 
that the PA of the Be disk is ~ 115° (PA of polarization is perpendicular to the scattering 
plane defined by the light source, the scatterer, and the observer) and the rotational axis of the 
Be star is at PA ~ 25°. As noted in §1, radio jets have been found in LS I +61° 303. Table 2 
shows PA's of the radio jets in literature. The directions of the jets seem to be biased around 



Aand the central wavelength of the Ha emission band to be 6563 A. There is no significant line 




(8) 



(9) 
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PA~120° during quiescent phases. If we can assume that the radio jets are perpendicular to 
the plane of the accretion disk around the compact object, the PA of the semi-major axis of 
the accretion disk on the projected sky would be ~30°. No X-ray pulsation (e.g., Goldoni 
& Mereghetti 1995) suggests that the compact star is not strongly-magnetized and such an 
orthogonal disk-jet geometry can be assumed. We also consider that the plane of the accretion 
disk represents the orbital plane of the binary system. Thus, we propose a geometrical model 
in which the rotational axis of the Be star is almost perpendicular to the orbital plane as shown 
in Figure 6. Considering an evolution of a typical binary system whose member stars were 
born in a same molecular cloud, our geometrical model seems to be quite peculiar. However, a 
similar geometry has been reported in another Be/X-ray binary system, PSR B 1259-63 (Wex 
et al. 1998 and Wang et al. 2004). They suggested that the Be equatorial disk is considerably 
tilted with respect to the orbital plane, and Chernyakova et al. (2006) suggested that the tilt 
angle is about 70°. These studies may support our picture for LS 1 +61° 303. We will discuss 
possible origin of this peculiar topology in §4.4. 

4-2. Inclination of Be equatorial disk 

Considering a light scattering model derived by Cassinelli, Nordsieck & Murison (1987), 
we can estimate the inclination angle i' of the rotational axis of the Be disk against the line of 
sight. The expected degree of the polarization is given by 
3crTsin^z' 



P= - 



/ / {l-3fi')D{r)n,{r,fi)dfidr (10) 

16 JrJfi 

where P is a degree of polarization, ctt is the Thomson scattering cross-section, /i is cos 6 and 
^ is a polar angle against the rotational axis of the Be star. D{r) = (1 - Rl/r'^f'^ is the finite 
disk depolarization factor. D{r) rises from to 1 as r changes from i?* (stellar surface) to +oo. 
rie is an electron number density of the ionized gas. We assume that the ionized gases exist in 
a disk form as 

f n,{r) = ^^~^}°''^ X — (90° - g,,, < g < 90° + 
ne{r,fi) = < Airv{r)r^ nin (11) 

[ (^ < 90° - Ohoa and 90° + Ohoa < 0) 

where is a wind mass loss rate, fihoa = cos 9hoa, dhoa is a half opening angle of the equatorial 
disk, mu is a mass of hydrogen atom, and v{r) = vq + {voo - vq){1 - R^/r"^) is a wind velocity 
of Be star (Castor & Lamers 1979). We set the initial velocity to vq ~ 5km s~^ and the 
terminal velocity to ~ 200 km s~^ (Waters et al. 1988). Adopting P ~1.5%, the intrinsic 
polarization at ~ 4000 Awhere the effect of the unpolarized bound-free emission of hydrogen 
is nearly neghgible, Ohoa ~ 7-15° (Waters et al. 1988; Marti' & Paredes 1995), and ~ 
(0.4-4) X lO-^Mo/yr (Marti & Paredes 1995), we can get i' ~ 20-50° from Equations (10) and 
(11). Ultraviolet spectroscopy indicates that the rotational velocity of the Be star is V sin i' = 
360 ± 25 km s~^ (Hutchings & Crampton 1981). By assuming that the maximum rotational 
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velocity of Be stars is Vmax = 630 x 0.9 = 570km s~^ (Hutchings et al. 1979), Massi et al. 
(2001) suggested the lower limit for the inclination of the equatorial disk as ~38°. Combining 
our result, we suggest that the inclination angle of the Be disk is i' ^ 38-50°. 

4.3. Geometrical model of LS I +61° 303 

Orbital parameters of LS I +61° 303 have been derived from optical spectroscopy as 
orbital eccentricity e = 0.72 ± 0.15 and a sin i = 8.2 ± 2.9i?0 where a is the length of major 
axis of the orbit and i is the inclination angle between the orbital axis and the line of sight 
(Casares et al. 2005). From e and a sin i, we can estimate the distance of the compact star 
from the Be star at their periastron passage as dpcri = a (1 - e). Considering the range of ± 
la error in both e and asini, we calculate rfperi for i = 0° - 90° as shown in Figure 7. This 
figure suggests that the separation at the periastron passage is less than ~ 3Oi?0 in most cases 
for i ^ 10°. If the outer radius of the Be equatorial disk is ~ 3/2* = 30Rq (Waters et al. 
1988), the compact star is likely to get across the Be equatorial disk, nearly perpendicularly 
as discussed in §4.1 For larger i (^ 30°), the compact star would pass through the inside of 
the Be stellar surface. In this case, a drastic dip of high energy emission would be observed 
around the periastron passage ($ ~ 0.23). However, such a phenomenon has not been reported 
(e.g. Paredes et al. 1999; Harrison et al. 2000; Massi 2004b), and we can reject the case of 
c^peri ^ R* and suggest i ^ 30°. This condition is comparable to the condition ( ? ^ 28 ± 11°), 
derived under the assumption that the compact star is not eclipsed by the Be star (radius ~ 
10 -R0) even at the conjunction ($ = 0.16, Casares et al. 2005). When the compact star gets 
across the Be equatorial disk, the accretion disk around the compact star may obtain angular 
momentum which is different from that of the spherical wind of the Be star. Using Bondi-Hoyle 
accretion modeP (Bondi & Hoyle 1944) together with mass-loss model in Waters et al. (1988) 
(woo ~ 2000km s~^, and ~ 1Q~^Mq yr~^ for the spherical wind, and ~ 200km s~^ and 
Mw ~ 1O~^M0 yr~^ for the equatorial wind), and a disk geometry model same as §4.1, we 
estimate the mass accretion rate on to the compact star. The accretion rate during the passage 
through the equatorial disk at ~ 3/2* (A/^^disk) is ~ 8 x 1O~^M0 yr~^ and the spherical wind at 
the same passage (M^^sw) is ~ 7 x 1O~^^M0 yr"^ Thus, M^^disk is about 6 order higher than 
Mw,sw The derived value of M^,disk is about 40 times higher than the Eddington accretion rate 
of a typical neutron star, and it is unlikely that such a high rate accretion actually occurs in the 
case of a neutron star primary. But, the large difference in Mu, would lead to sudden changes 
of accretion properties and the resultant mass ejection. Such an impulsive change gives a hint 
on the formation and the precession of the radio jet. However, the launch of the precessing 
jet occurred at $ ~ 0.7, largely apart from the periastron $ ~ 0.23. It has been known that 
both the radio and X-ray emission modulate with the orbital phase (Harrison et al. 2000 and 
references therein) and that there is a significant offset between their peaks (A$ ~ 0.4-0.5). The 

^ Assuming Kepler rotation, we found the rotational velocity of Be equatorial disk (250kni s~^). 
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radio peaks are distributed over a wide phase interval of about 0.45-0.95 (Paredes et al. 1990). 
Harrison et al. (2000) suggests that X-ray emission is produced by inverse Compton scattering 
of stellar photons by relativistic electrons produced at the shock interaction of the relativistic 
between jet and Be star wind near the periastron, and the synchrotron radio emission would 
then arise from the expansion of the plasma around the apoastron passage. Along with this 
scenario, the launch of the precessing jet at $ ~ 0.7 does not seem to be directly connected with 
the activities inferred from orbital modulations of X-ray and radio fluxes. More speculatively, 
we suggest that the time lag between the periastron passage and the jet launch represent the 
timescale of the mass accretion through the accretion disk. To examine the mechanism of the 
jet formation we should gain more observational information, e.g., more precise X-ray and 7-ray 
monitoring around $ ~ 0.7. 

4.4- On the origin of peculiar topology 

Mirabel, Rodrigues & Liu (2004) suggested that the past supernova explosion in LS 
I +61° 303 (which produced the compact object) should be asymmetric to kick the compact 
object up to the highly eccentric (e ~ 0.7) orbit. The asymmetric explosion would also give 
a proper motion to the binary, fi ~ 1.55 mas yr~^ (27 km s~^) and its PA ~ 141°, escaping 
from its possible birth place IC 1805 (Lestrade et al. 1999; Mirabel, Rodrigues & Liu 2004). 
X-ray observations of several young pulsars have provided evidence for approximate alignment 
between the pulsar proper motion and its spin axis (Lai et al. 2001; Romani & Ng 2003). 
This may lead the spin axis of PA ~ 141° for the compact object of LS I +61° 303, if it is 
a neutron star. However, it is unclear that we can adopt the spin-kick alignment mechanism 
to this binary because the proper motion is much less than the typical kick velocity of several 
hundred km s~^. The asymmetric explosion might have had effects selectively on the binary 
orbit only. 

5. Summary 

We performed optical spectropolarimetry of the microquasar LS I +61° 303. We precisely 
estimated the ISP component from the Ha line emission and derived the intrinsic polarization. 
Our main results can be summarized as follows: (1) The wavelength dependence of the intrinsic 
polarization is explained by Thomson scattering in the equatorial disk of Be star. (2) We found 
that the PA of the semi-major axis of the Be disk on the projected sky is almost parallel to the 
PA of the radio jet observed in quiescent phase. If we assume that the radio jet is perpendicular 
to the accretion disk of compact star, the Be disk is approximately orthogonal to that of the 
accretion disk. (3) The inclination of the Be disk is i' ~ 38-50°. (4) The compact star is likely 
to get across the Be disk around the periastron passage suggesting a drastic change of mass 
accretion rate (Mw,disk/^«;,sw ~ 10^). Speculatively, we suggest that this might be related with 
the precession of the radio jet axis. LS I +61° 303 is one of the the ideal laboratories to research 
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the peculiar activity of microquasars. Future phase-resolved multi-wavelength monitoring is 
required to examine the jet launching. 
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Fig. 1. Polarization spectra of LS I +61° 303 on 2005 Oct 27 and Nov 8. Upper panel (a) shows degree 
of polarization Pobs, and the lower one (b) shows its position angle ^obs- The polarimetric data are binned 
to a constant photon noise of 0.02%, and the observational error(lcr) in each bin is shown by a vertical 
error bar. 
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Table 1. Observation log 
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EW of Ha 


^obs 


^obs 


rxi 


nd 
%a 








(s) 


(A) 


(%) 


/o \ 

( ) 


(%) 


( ) 


2005 Jan 21.5 


2453392.0 


0.37 


300x36 


7.98±0.09 


1.22±0.03 


138.8±0.7 


2.21±0.93 


129.3±18.6 


2005 Jan 28.5 


2453399.0 


0.63 


300x44 


9.89±0.15 


1.31±0.03 


137.8±0.6 


1.85±0.51 


127.6±14.7 


2005 Oct 27.6 


2453671.1 


0.90 


300x64 


9.61±0.09 


1.20±0.02 


137.2±0.5 


2.63±0.36 


125.6±7.00 


2005 Oct 30.7 


2453674.2 


0.94 


300x68 


9.15±0.11 


1 Of\-l-ri no 
i.zDitU.Uz 


1 QQ 0-\-ri K 

ioo./±U.o 


/.oo±U.4 / 


izo.oiniU.i 


2005 Nov 6.7 


2453681.2 


0.28 


300x40 


9.30±0.11 


1.30±0.03 


137.1±0.6 


2.53±1.01 


130.9±20.1 


2005 Nov 7.7 


2453682.2 


0.32 


300x68 


lO.lOiO.ll 


1.31±0.04 


137.9±0.8 


1.74±0.57 


128.0±17.0 


2005 Nov 8.6 


2453683.1 


0.35 


300x56 


lO.lliO.ll 


1.27±0.02 


137.8±0.4 


1.67±0.37 


128.3±11.5 


2005 Nov 10.6 


2453685.1 


0.43 


300x60 


9.49±0.13 


1.25±0.02 


135.9±0.4 


2.44±0.92 


117.8±20.6 


' Orbital phase calculated for to = 


JD 2 443 366.775 and P= 


=26.4960days 


(Gregory 


2002). Total 


exposure 



time is expressed as the integrated time per one frame, multiplied by the number of frames. degree of observed 
polarization and it's position angle in the synthetic V-band (Bcssel 1990; Kawabata et al. 1999). Polarization 
of the Ha line emission component. See §3.3. 



Table 2. Observation of radio jet in the literature 



Epoch 


Array 




State 


PA 

n 


1987 Sep 25" 
1987 Oct 1" 


EVN 


0.56 
0.79 


quiescent 
burst (decaying) 




1990 Jun 6'' 


EVN+VLA 


0.71 


burst (decaying) 


~135,~30 


1992 Jun 8^ 


Global VLBA 


0.39 


quiescent (minioutburst) 


-160 


1993 Sep 9'' 
1993 Sep 13^* 


Global VLBA 


0.66 
0.81 


burst (variable) 
quiescent 


-120 


1994 Jun 7^= 


EVN 


0.89 


quiescent 


-120 


1999 Sep 16-17^ 


HALCA+Global VLBA 


0.65 


burst (variable) 




2001 Apr 22-239 


MERLIN 


0.68 
0.71 


quiescent 
preburst 


-124 
-67 



"Taylor et al. 1992; ''Massi et al. 1993; ^Peracaula et al. 1998; ''Paredes et al. 1998; ''Massi et al. 2001; 
•^Taylor et al. 2000; ^Massi et al. 2004a; ''Orbital phase is calculated in the same method of Table 1. 
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Fig. 2. Enlarged plot of Fig.l around the Ha line. From top to bottom, we plot (a) relative flux (flux at 
6500 A=: 1), (b) degree of polarization and (c) its position angle. In (a) airmass effects are not calibrated, 
but the data are corrected for the instrumental response function. We can see a clear change of polarization 
across Ha emission line, which suggests an existence of a polarization component intrinsic to the target 
star. 
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Fig. 3. Equivalent width of Ha emission line as a function of the phase of mass- loss determined by radio 
flux (Gregory ct al. 1999). Crosses are from Zamanov et al. (1999), inclined crosses are Parcdes et al. 
(1994), and circle denote our results. Our results are commonly small compared with past observations. 
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'^max (angstrom) 



Fig. 4. K vs Amax derived for our nightly-averaged data (open square) are plotted together with samples 
of typical ISP (filled square). The sample stars are from Whittet et al. (1992). The solid line is K = 
1.66Aniax (/iin) -|- 0.09, which well represents the correlation between K and Amax of ISP (Whittet et al. 
1992). The K values of our data are commonly apart from the line (except for 2 data points, which would 
be affected by noise), which indicates that the observed polarization cannot be explained by only ISP. 
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Fig. 5. Spectra of intrinsic polarization, (a) degree of the intrinsic polarization and (b) its position angle. 
The data are binned with the same manner in Figure 1. The wavelength dependence of polarization is 
quite similar to that of a typical Be star, e.g., C Tau (Wood ct al. 1997). In a Be star, intrinsic polarization 
is originated by wavelength-independent Thomson scattering in the equatorial disk. The gradual decrease 
of polarization with wavelength from ^4000 to ~8000 A and the polarization jump beyond 8000 A are 
explained by unpolarized bound-free emission from hydrogen plasma. 
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Fig. 6. Schematic representation of circumstellar geometry model of LS I +61° 303. The rotational axis 
of the Be star is nearly perpendicular to that of the accretion disk around the compact star. The compact 
star is likely to get across the equatorial disk of the Be star around its pcriastron passage (see figure 7) . 
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Fig. 7. Estimated binary separation at the pcriastron passage (dpori) as a function of orbital inclination 
i. For the calculation we used orbital parameters derived by Casares et al. (2005) with typical stellar 
parameters of a BOV star. The range of probable dpcri values lies between the two dashed curves. The 
compact star is likely to get across the equatorial disk of the Be star. 
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